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Primordial Lepton Family Asymmetries in Seesaw Model 
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In leptogenesis scenario, the decays of heavy Majorana neutrinos generate lepton family 
asymmetries, Y e , Y M and Y T . They are sensitive to CP violating phases in seesaw models. The 
time evolution of the lepton family asymmetries are derived by solving Boltzmann equations. 
By taking a minimal seesaw model, we show how each family asymmetry varies with a CP 
violating phase. For instance, we find the case that the lepton asymmetry is dominated by 
Y,j, or Y T depending on the choice of the CP violating phase. We also find the case that the 
signs of lepton family asymmetries Y^ and Y T are opposite each other. Their absolute values 
' can be larger than the total lepton asymmetry and the baryon asymmetry may result from 

| the cancellation of the lepton family asymmetries. 



< 



■ §1. Introduction 

*sO ' 

\ From the recent astronomical observations, the precise values for cosmological 

qq ' baryon to photon ratio (f? = — ) are obtained, 

o' 

g; 7 ? =(6.5±°;|)xl0- 10 E (6.0±Ji)xlO- 10 0. (1-1) 

Leptogenesis^ is a very attractive scenario for the matter and the anti-matter asym- 
Ph. metry of our universe.™ Because the seesaw mechanism^ which underlies the sce- 

Q_i| nario is able to give a natural explanation for tiny neutrino masses, it is important 

to study the scenario from both cosmological and particle physics point of views. 
In the early literatures EU3 the lepton number production and the evolution in the 
expanding universe were studied. In the last few years, much attention has been paid 
^ to CP violation of seesaw models. The relation between CP violation of leptogenesis 

and CP violation of neutrino oscillation was shownP Because there are many CP 
violating phases in seesaw models, it is important to identify CP violating observ- 
ables as much as possible. Such observables include the low energy CP asymmetry 
of neutrino oscillations. In the present paper, we study the lepton family density 
to entropy density ratios Y e , and Y T in seesaw models. Hereafter, we call them 
the lepton family asymmetries. We shall show that they are sensitive to the CP 
violating phases. By studying the lepton family asymmetries, we can trace more 
detailed history of the leptogenesis scenario. 

So far, lack of our knowledge for the neutrino Yukawa couplings in seesaw models 
has resulted in numerous speculations on their possible forms. Our result shows that 
different forms may lead to different lepton family asymmetries. Thus, the lepton 
family asymmetries are useful for constraining the forms of the Yukawa couplings. 



E-mail: endoh@theo.phys.sci.hiroshima-u.ac.jp 
E-mail:morozumi@theo. phys.sci.hiroshima-u.ac.jp 
E-mail: xiongzh@theo.phys.sci.hiroshima-u.ac.jp 



typeset using ^TpTeX.cIs (Ver.0.89) 



2 



T. Endoh, T. Morozumi and Z. Xiong 



In this paper, we give a detailed calculation for the lepton family asymmetries 
by solving the Boltzmann equations. In the Boltzmann equations of the previous 
works P ™ it is implicitly assumed that the chemical potentials are the same for 
all the three lepton doublets and the effects of the chemical potentials of the other 
particles are neglected. Contrastingly, we shall consider all the chemical potentials 
for the standard model (SM) particlesP We do not assume that the three lepton 
doublets have the same chemical potential a priori. We assign the independent 
chemical potentials for each generation of lepton doublets.^ In this way, we can 
treat the evolution of the lepton family asymmetries in the Boltzmann equation. 

For numerical analysis, we extend the previous studjCU on the minimal seesaw 
model. In the model, there are three CP violating phases. Considering the lep- 
ton family asymmetries, we find that they are sensitive to one of the CP violating 
phases which is different from the one related to the total lepton asymmetry Y^. 
For instance, we find the case that the signs of the lepton family asymmetries Y^ 
and Y T are opposite each other and the tiny lepton asymmetry Yl results from the 
cancellation among the much larger family asymmetries. 

The paper is organized as follows. In Sec. |2J the Yukawa sector for seesaw models 
is given and the minimal seesaw model is introduced. In Sec. |3 The Boltzmann 
equations for the lepton family asymmetries are derived. Sphaleron process and the 
other equilibrium processes are discussed and the relations among chemical potentials 
are studied. Using the minimal seesaw model and the relations of the chemical 
potentials, we solve the Boltzmann equations and present the numerical results in 
Sec. 0J Finally, in Sec. we give our conclusions. The derivations of various cross 
sections used in our calculation are shown in Appendices. 

§2. Model 

We first discuss the leptonic sector of seesaw model. The Yukawa terms and the 
mass terms for seesaw model are given by, 

C m = -yl%N Rk 4> - y\Lll Rl 4> ~ \n^M Nr N r + h.c, (2-1) 

where i = 1,2,3 and k = 1, 2, • • • N. Li are SU(2) lepton doublet fields, Nn k are 
the heavy Majorana right-handed neutrinos and l Ri are the right-handed charged 
leptons. M]\r R is NxN Majorana mass matrix of the right-handed neutrinos and is 
a diagonal matrix, i.e., M^ R = diag.(Mi, M2, ■ ■ ■ , Mn). yi are the Yukawa terms 
for charged leptons. We can take the basis in which both and y\ are real and 
diagonal without loss of generality. In this basis, flavor violating processes occur 
through off-diagonal elements of the 3xN Yukawa matrix y v . In the broken phase, 
Higgs field has the vacuum expectation value v = 246 [GeV], and Dirac mass term 
is generated as mp = -^mVv 

The minimal seesaw model which is compatible with the present neutrino os- 
cillation data includes two heavy Majorana neutrinos. For numerical analysis, we 
focus on the minimal model. We summarize the results in the previous worli^ 
which are relevant to this paper. There are 11 parameters in neutrino sectors in 
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the model. They are two heavy Majorana masses, M\,M 2 and 9 parameters in 
rriD. As input parameters, we choose a heavy Majorana mass M±, the ratio of the 
two heavy Majorana masses, R = and their total decay widths rj^(k = 1,2). 
From neutrino oscillation experiments, two mass squared differences and two mixing 
angles are obtained. Then 8 of 11 parameters can be determined. The remaining 
three parameters are left undetermined. They are related to \U e z\, CP violation of 
neutrino oscillation and neutrinoless double beta decays, which will be measured in 
the future experiments. The formulae which relates the input parameters with mo 
are derived using the bi-unitary parameterization,^ 



m D = U L mV R , 



(2-2) 
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TTi.2 and m 3 are real and positive and m 2 <m 3 . Notice that there are three CP 
violating phases Sl,Jl and jr. 

The four parameters in m and Vr can be written in terms of the heavy Majorana 
masses M k , their decay widths and the light neutrino masses rii(i = 1,2,3), 
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where we define the auxiliary quantities as, 

{m ] D m D ) kk 
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The above relations follow from the eigenvalue equation for the light neutrino mass 
squared n 2 , 

.2\ 
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m 



eff 



-m D M^m T D . (2-7) 



We can easily see that one of the eigenvalues {n\) is zero. It was shown that MNS 
matrix can be written as a product of matrices as follows, 



Umns = U L K R , (2-8) 



where Kr is given by, 



\ / 1 

K R = I cos6» sinfle-^ e ia ). (2- 



smtfe 



cos 6 / \ 



6, 4> and a in Kr are functions of the parameters in m and Vr which are already 
determined by Q2-5j) . Explicitly, Kr is obtained from the diagonalization of Tn e ff as, 

Ul NS m ef fU* MNS = -K R (mV R M^V^m T ) K* R = diag.(0, n 2 , n 3 ), (2-10) 

To summarize this section, we are able to determine the part of MNS matrix, i.e., 
Kr using the parameters Mi, R, xi, x%, n 2 and n%. At this stage, Ul is still left 
undetermined. How to constrain Ul will be discussed in Section 0J 

§3. Boltzmann equation 

In this section, we shall give the derivation of the Boltzmann equations for 
lepton family asymmetries. The primordial lepton numbers are generated by out- 
of-equilibrium decays of the heavy Majorana neutrinos. The decay occurs at the 
temperature higher than the electroweak phase transition temperature (Tew)j be- 
fore sphaleron process being frozen. While at extremely high temperature T > 
Tsph ~ 10 12 [GeV], the sphaleron process is not so active and the lepton number 
can not be converted into the baryon number. Then, the primordial leptogenesis 
must occur at the temperature above Tew and below T sp h so that the baryogenesis 
from leptogenesis works. Above Tew> the electroweak symmetry is recovered and 
the Higgs vacuum is in the symmetric phase (v = 0). We focus on the Boltzmann 
equations which are valid in the symmetric phase. In the phase, the lepton family 
number violating processes involve the Yukawa coupling of the neutrinos y v . Using 
the Boltzmann equation, we can trace the evolution of the lepton family number 
asymmetries from the mass scales of heavy Majorana neutrinos down to Tew- I n 
order to trace the evolution into the lower temperature regime T < Tew> the Boltz- 
mann equations in the broken phase must be used. The derivation of the equations 
in the broken phase is beyond the scope of the present paper. 

We first show the Boltzmann equations for the heavy Majorana neutrino number 
densities (k = 1, 2, ■ • ■ N) and lepton family number densities (i = 1, 2, 3), 



dn k 

dt 
dn L 
dt 



+ 3H(z)n k = C nk , 

+ 3H(z)n Ll =C k , (3-1) 
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where Mi is the mass of the lightest heavy Majorana neutrino and z = -ijr. H(z) = 

H(l)z~ 2 is the Hubble parameter with H(l) = \J ^l®* mp denotes the Planck 
mass scale and ~ 106.75 counts the total number of light particles degrees of 
freedom at Tew < T < T sp h- The right-hand side of the Boltzmann equations C nk 
and C\ i are source terms and they denote the particle number changes per unit 
spacetime volume. The source terms are obtained by computing the decays, the 
inverse decays, and the scattering processes. 

In deriving the Boltzmann equations, we use Maxwell-Boltzmann distributions 
for light particles. Explicitly, the phase space density of a light particle X with 
energy Ex and chemical potential fix is given by, 



fx = exp 



Ex — Hx 



(3-2) 



The distribution is a good approximation in the absence of Bose condensation and 
Fermi degeneracy. The signs of the chemical potentials for particle and anti-particle 
are opposite each other. The up component in SU (2) doublets has the same chemi- 
cal potential as that of the down component in the symmetric phase of electroweak 
symmetry. For the heavy Majorana neutrino Njt k , we assume that the distribution 
is proportional to the equilibrium distribution. Under this assumption, the normal- 
ization constant is given by the ratio of the non-equilibrium number density and 
the equilibrium density n^ q , 



f N (s z = h = Ms, = -b = ^e~ E -/ T , (3-3) 

— ■ Zi TX-i 



k 



eq 



where s z denotes the spin of the heavy Majorana neutrinos, = M^/Mf. K\ and 
K2 are modified Bessel functions. 

3.1. Calculation of C nk and 

The processes contributing to the source terms can be divided into two parts. 
One includes the decays and inverse decays of the heavy Majorana neutrinos, and the 
other includes two-body scattering. The contribution from decay and inverse decays 
are denoted as C D and the contribution from scattering processes are denoted as 
C R . Then the total source terms are given by their sum, C nk = C® k + C„ k for the 
heavy Majorana neutrinos densities and = + for the lepton family number 
densities respectively. 

Below, in the results of the two-body scattering a + b — > c + d, we define the 
reaction rate 7 as, 



Mf 1 

64vr 4 i 

(m a +m b y/Ml 



00 

J dsa(s)ViKi(zV§), (3-5) 
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where a(s) is the reduced cross section and s = (p a + Pb) 2 /Mf. The reduced cross 
section a(s) is related to the usual cross section a(s) by a(s) = j- + pb y2 [(Pa • Pb) 2 — 
m 2 .ml]a(s). All the reduced cross sections a are collected in Appendix B. 



(a) 



N Rk 



(b) 



I- 



'J 



N Rk 



+ 



(d) 



Fig. 1. Feynman diagrams of the heavy Majorana neutrino decay: (a) Tree level diagram. One 
loop diagrams: (b) vertex type and (c) (d) self-energy type diagrams. 



Firstly, we calculate the source term 
C n . k in (|3-lj) . The decays and the inverse 
decays of the heavy Majorana neutri- 
nos contribute to the source term. The 
processes involved in the calculation are 
listed in Table [IJ Using the Lagrangian 
(|2-1|) . the definitions in Appendix A, 
and the phase space densities in (|3-2j) 
and ()3-3|) . we calculate the Feynman di- 
agrams shown in Fig. Wi a ) and obtain, 

n D ST^ ki ( n k 



Since the Yukawa coupling of top 
quark is much larger than those of the 
other SM particles, the heavy Majorana 
neutrino creation and annihilation pro- 
cesses in Fig. 12 give important contri- 
butions to the heavy Majorana neutrino 
number densities.® The processes are 
listed in Table [HJ Their contributions 
to the source term C nk are given by, 

< = C#* + Cg)*, (3-8) 



Table I. The decays and inverse decays of 
the heavy neutrinos contributing to lepton 
and heavy Majorana neutrino number den- 
sities. 



Number change 


Processes 




AL = 1, AN = 1 


N Rk - 
N Rk <-> Vi(j> 
N Rk <-> Vi<tP* 


Id 


source term 








««„ (a) 5 (b) 



Fig. 2. Processes contributing to both the 
heavy Majorana neutrino densities and lep- 
ton number asymmetries with Higgs ex- 
changes via (a) s channel and (b) t channel. 



where 7^ is related to the partial decay width r k i{N Rk —>■ l i (p + ), 
7% = 4nl«r ki §¥^B r kl = -L \{y u ) lk \ 2 M k . 
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Table II. The processes with Higgs exchange contributing to the lepton family number densities 
and heavy Majorana neutrino number densities. 



Number change 


Processes via t channel 




Processes via s channel 




AL = 1, AN = 1 


N Rk t R <-> irl L 
N Rk b L «-» irt R 
l+b L «-> N Rk t R 
lft R » N Rk b L 
N Rk t L «-> tij^i 
t R Vi <-> N Rk t L 
N Rk t R <-> Tl^i 
iL^i <-> N Rk t R 




&z7fl «-» N Rk l~ 
thin <-» N Rk Vi 
N Rk Vi <-> t^ta 




source term 











where C^^ and Cn} R 



arc, 



^ (cosh ^+ cosh ^ 



, ^- + ^ , Mi- 
cosh — — cosh ■ " ' 



/7(2)« 



—go — cosn 



T 



(3-9) 
(3-10) 



The reaction rates 7$* „ and 7^ are related to the reduced cross sections <7?T„ and 



<T^* f . Their explicit form are given in (j±3- 15|) and (jij- 161) . The corresponding processes 
are shown in Fig. HJa) and Fig. EJb), respectively. 

Next let us study the source terms for lepton family numbers, C\. in Q3TJI . The 
decays and the inverse decays contribution to Cf is given by, 



riD \^ hi 
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(3-11) 



where the lepton family CP asymmetry e\ is defined as 



r(N E 



1: 



r(N Rk - i+<f>- 
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(3-12) 



The CP asymmetries are generated by the interference of the tree diagram and one- 
loop diagrams shown in Fig. ^ 
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where Xk'k = M%,/M% and 
I(x) = 



H \-(l + x)ln 



-3--1/2 

-2x 3 / 2 



1 — X 1 + X 

for x S> 1, 
/or s<l. 



(3-14) 



The diagram Fig. ^d) gives the contributions to the second terms of the lepton 
family CP asymmetry (|3-13|) , The diagram does not contribute to the total lepton 
asymmetry while it does contribute to the lepton family CP asymmetries. 

The contributions to the source term from two particle scatterings can be 
divided into four parts, 



(x)R 



(3-15) 



x=i 



The first two terms in (|3-15|) are obtained by calculating the diagrams in Fig. 2. 
They are given as follows, 



CH )R [Fig. Etb)] = £ 7ft fe (smh ^ - sinh 



T 



+ smh lL — hsmh ' 



cf*[Fig. 2(a)] = £t£ 



T 



n 



sinh — =^ + sinh 



r 



(3-16) 
(3-17) 



The processes listed in Table III also contribute to the C^. They are denoted by 
Cj and and are given by, 



Table III. \AL\ — 2 and AL = lepton family number changing scattering processes 



Number change 


Processes via s (u) channel 




Process via t (u) channel 




\AL = 2\, AN = 


If ^ - If^ 

Vi4>°* <-> i/ 3 >° 




U *0 U * <-> I/jl/j 

2f 2+ «-> 0±0± 

0°0° «-» 


7^,2 


F^ u * ^ l-<f> + 
1+$- <_> Vi< f 

lf<t>- <-» 
77,/* <-> 27 0+ 




27^ <-> 0~0 U * 
2<Fj «-> 0+0° 

IJVi <-> 0~0°* 
2+77, 0+0° 




AL = 0, ALt ^ 


if 0+ - if 0+ 
VjcfP* <-> 2+0" 
Vi(f>°* <-> Vj4>°* 

v j( j> <-> l~4> + 

77 s 0°* <-> 270+ 
270+ ^0° 




0+0" « If 2+ 

ViVj «-> 0°0°* 

VjVi <-> 0°0°* 
Uilf <-> 0+/* 
77727 «-> 0°0~ 
777/7 <-> 0°0~ 
0+0°* 


7iV,£2 


source term 
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(3-19) 




\ 



0+ 



(a) 



(b) 



Fig. 3 



The lepton family number violating 
scatterings with Nn k exchanges via (a) s 
channel and (b) u channel. 



The processes in Table III are classified into the total lepton number changing scat- 
terings (|AL| = 2) and the total lepton number conserving scatterings (AL = 0). 
In the former, the lepton family number 
Li changes by one unit or by two units. 
The latter corresponds to \ALi\ = 1. 
The total lepton number conserving and 
the lepton family number violating pro- 
cesses must be taken into account when 
we study each lepton family number 
evolution. In Fig. 3, the typical 
Feynman diagrams which correspond to 
AL = ALi = — 1 are shown. Their 
reaction rates are included in "Yn,s and 
7iV,t2 of ()3-18j) and (|3-19|) . Finally we also note that in (|3-18j) . the term e\ comes from 
the on-shell contribution of s-channel heavy Majorana particle exchanged diagrams. 

3.2. Sphaleron and chemical potential relations 

In the previous subsection, we write the source terms with the reaction rates 7 
and the chemical potentials /x. In this subsection, we derive the relation between the 
lepton family asymmetries Yn and chemical potentials by considering the various 
equilibrium conditions. We adapt the method developed by Harvey and Turner. 

At first, let us consider the effective action of sphaleron for all left-handed 
fermions, 

O b+ l = H(Q l QiQ i L i ), (3-20) 

i 

where Qi,Li are quark and lepton doublets respectively. The sphaleron transition 
rate at high temperature, in the symmetric phase, is estimated to be 121 

2xl0 2 a^T, (3-21) 



sph 
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where a\y is weak coupling constant at Grand Unified scale and is taken to be 
jq. From the thermal equilibrium condition, i.e., -T sp h > H , we obtain the upper 
bound of the temperature T sp h = 1.4 x 10 12 [GeV]. Above the temperature, sphaleron 
process is not active. At low energy, below Tew — 10 2 [GeV], the sphaleron process 
is frozen out again, then the active temperature range of sphaleron is, 

Tew <T< T sph . (3-22) 

In the temperature range (|3-22j) . the sphaleron equilibrium condition leads to the 
chemical potential relation, 

N f 

£ (»u iL + 2f*d lL + IHu) = 0. (3-23) 
i=l 

where Nf = 3. We assume that gauge and Yukawa interactions are in equilibrium. 
From charged current interaction processes, we obtain the following chemical poten- 
tial relations, 

Pw- = Vd lL - Vu iL = fi r ~ Vu iL = (if - + /tyo, (3-24) 

where i (i = 1,2,3) denote indices for the generations. By taking the weak basis 
that the Yukawa term for up type quarks is flavor diagonal and that for down type 
quarks has flavor off-diagonal elements, we obtain the chemical potential relations 
from the equilibrium condition of the Yukawa interactions, 

H° = Vum - Vu lL = lid iL ~ Vd jR = Mr r ~ h~ ' ( 3 ' 25 ) 

where we note that the chemical potentials for the left-handed down quarks {dn) 
and those of the right-handed down quarks {djRj satisfy the flavor mixed relations. 
Therefore the chemical potentials for the right-handed down quarks and the left- 
handed down quarks are flavor independent. Both of them can be written by a 
single chemical potential as, 

Hd iR = Vd R , Vd iL = Wl' (3-26) 

Using l|3-24|) . (|3-25|l and (|3-26|) . we can also show that the chemical potentials for up 
type quarks are flavor independent, 

Mujij = M«h' ^ u iL = Mui- (3-27) 

Next, we can relate the chemical potentials of SU(2) doublets. This follows from the 
chemical potential for W boson vanishes in the symmetric phased From ()3-24|) . 
the up and down components of SU(2) doublets have the same chemical potential. 
Note that the chemical potentials for leptons can be flavor dependent. 



iL iR 



(3-28) 
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We also take into account of the charge neutrality condition. The condition can be 
written as, 



T 2 



^2 Vu iL - M - (2N f + N H + 2) fi w - + (2N f + N H )H4P 



i=i 



0, 



(3-29) 



where fj,= Hv iL an d Njj is the number of Higgs doublet and is taken to be 1. With 

i 

(|3~23|) . (|3~26|l . (|3~27j) . (j3~29|) and fiw = 0, we obtain the following relations, 

9/i Mi + n = 

=fJ>~ 7^o. (3-30) 

Then the chemical potentials for quarks and charged leptons can be written by the 
single chemical potential of [i = ^ [i ViL of neutrinos as follows: 



(3-31) 



/'• 



We can also write the baryon number density and lepton number densities with the 
chemical potentials as, 



til 



n B 



T\ . T 2 ( 4 \ 



T 2 f 4 \ 17 2 

i=l v 7 



(3-32) 



T_ 
~6 



2 N f /2 \ 2 

i=l \ / 



By using the equations above, we finally obtain the relations between the chem- 
ical potentials for the lepton doublets (jjl^- ) and lepton family asymmetries, 





1 


s 


~ 2 


TIL 


17 


s 


~ 42 



s 51 



(3-33) 



where s is entropy density given by s = ^-g* (^) • What is important particularly 
from the flavor-dependent point of view is that family asymmetry Yjj t take a different 
value on the each generation. The conversion rate from lepton asymmetry to baryon 
asymmetry is given by the well known formulae,^ 

no 

Y B = ~Y L . (3-34) 
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We will use the relations (|3-31[1 and (|3-33[1 to express all the chemical potentials in 
terms of the lepton family asymmetries in the next subsection. 

3.3. Boltzmann Equations 

Now we can write the Boltzmann equations in a tractable form. Using the 
definition in (|3-33|) .y;v. = ^f- and Y^ q = the Boltzmann equations (|3-1|) can be 
rewritten as, 



dY Nk 
dz 

dY Ll 
dz 



z Y Nk 



1 



E b k D + + 27& 



aff(l) 



E 



Y N m 



(3-35) 



where the time t derivative in (|3-1|) is replaced by derivative on z using the relation 
in the radiation-dominated epoch, 



1 



2H 1 



H(l) 



45 mp 



(3-36) 



In the source terms C®' 11 , Cn k ,R derived in subsection 13. 11 we adapt the chemical 



potential relations in 



and (|3-33|) . and the approximation of sinh 



cosh ^ ~ 1. Then Aij and are given as follows, 



,4, 



16 E 



4 \ 



{Sij + 5 jf ) + _ W 7j * 1 + ) + ( — + — (4, + ) 7 ^ 2 



1 



51 16 



E 



1 - ^ \ ki 1 ki 



17 



ki 



17 



+ <&« + 



51 



(3-37) 
(3-38) 



§4. Numerical results 



In this section, we present the numerical results for the lepton family number 
asymmetries based on the minimal seesaw model described in Sec. El As we showed 
in the previous section, the Yukawa coupling y u = V ^ ,J) , and the heavy Majorana 
masses M%,M2 are needed for the numerical analysis of Boltzmann equation. We 
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parameterize mo in the bi-unitary form, i.e., mo = UimVn and write the parameters 
in m and Vr with M\,R = j^,xi,x 2 ,n 2 and 713. 

Next we show how the angles in Ul are constrained by discussing the presently 
available neutrino oscillation experimental measurements. The SK Collaboration 
showed that the created in the atmosphere oscillates into v T with almost max- 
imal mixing,^ sin 2 20 a tm ~ 1 and the neutrino mass squared difference Anv\ tm = 
(2 ~ 4) x 10" 3 [eV 2 ]. The second mass-squared difference and mixing angle are 
constrained by solar neutrino experiments. The SNO collaboration reported that 
the u e 's from the sun oscillate into the other active neutrinos™ The SK and 
the SNO collaboratiorP^ suggested that the MSW large-mixing-angle (LMA) so- 
lution is the most favorable solution to the solar-neutrino deficit problem, for which 
sin 2 26 so i = 0.7 ~ 0.9 and the combined analysis of the KamLANE^ and all the 
solar neutrino data gives Am 2 sol = (3 ~ 15) x 10~ 5 [eV 2 ]. For the third mixing angle, 
only the upper bound is obtained from the reactor neutrino experiments. CH002P^ 
found sin 2 29 rea < 0.1 for Z\m 2 tm ~ 3 x 10~ 3 [eV 2 ]. The current neutrino experimen- 
tal data indicates clearly that there is a hierarchy of neutrino mass-splitting. If the 
small mixing angles 9li3 and 9 are taken, the light neutrinos mixing matrix Umns 
defined in 1)2-8(1 . can be simplified, 



Umns — 



COS0L12 SU10L12 . . _ id) , 

- sin 9l\2 cos 9 L2 3 cos 9 L u cos 9 L2 3 sin 9 L2 3 

\ sin 9 L i2 sin 9 L2 3 - cos 9 L u sin 9 L2 3 cos 0^23 / 



P(a' 



(4-1) 



where <f)' = (ft + jl, a' = a — 7l/2 and P(a') = diag.(l, e ia ' , e~ ia '). A very similar 
form to the low energy MNS mixing matrix is obtained. In this case, the angles 
in matrix Ul defined in (|2-3() can be related directly to the corresponding neutrino 
mixing angles and can be determined by neutrino experiments. In this work, we take 
the natural hierarchical scenario and set the following neutrino masses in (|2-10|) for 
the corresponding measured neutrino mass differences, 

n 2 = ^Am\ ol = 7 x lO" 3 [eV], n 3 = ^ Am\ tm = 5 x 10~ 2 [eV]. (4-2) 
For the angles in Ul, we can take, 

7T 7T 

0L12 = 'sol = 7") 0L23 = 6atm = "T- (4-3) 
D 4 

The parameters x\ i2 defined in (|2-6(l are constrained by the conditionsp^ 

\x\ — x 2 \ < 113 — n 2 , n 3 + n 2 <xi+x 2 . (4-4) 

The primordial lepton number must be created when the sphaleron process is in 
equilibrium so that the conversion from the lepton numbers to the baryon number 
occurs effectively. In our numerical calculation, we fix the mass for the lightest heavy 
Majorana neutrino M\ as 2 x 10 11 [GeV], below T sp h ^ 10 12 [GeV] and we take the 
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ratio R = || to be 0.1. Finally, we consider the constraints from measurements of 
the cosmological baryon to photon ratio r] = n^/n 7 P^ Recall that temperature be- 
low Tewj sphaleron process is not active and thus, the total baryon number B is con- 
served. 

If the universe expands adiabatically, the 
total entropy S is also conserved. Then 
Yb is also conserved as, 

Y B (T EW ) = Y B (T = 3K) ~ If,. 



(4-5) 

Substituting (|4-5j) into (|3-34[) and com- 
bining the measurements in Q1-1JI . we in- 
fer the bounds on the primordial lepton 
asymmetry Yl = n^/s that we have to 
generate, 

Il(Tew) = ~n 

= -(1.15 - 2.13) x 10~ 10 , (4-6) 




Fig. 4. The evolution of the heavy Majorana 
Neutrino densities. 



at 90% C.L. 

We take the following values for the other 
parameters; 

sin6» L i 3 e" i<5L = 0, (x u x 2 ) = (0.0523, 0.0100) [eV]. (4-7) 

Concerning Higgs mass mjf, we change it from 150(GeV) to 800(GeV) and the results 
are not sensitive to the choice. So we show the figures for mu = 800 [GeVj. With 
them, only the CP violating phase jl is the parameter which remains to be fixed. 

To solve the Boltzmann equation, we start at T > Mi (z = 10~ 2 ) with the 
initial conditions, 



Ynh=Y%, Y Li = 0. (4-8) 

The typical results are displayed in Fig. to Fig- EI Fig. El shows the dependence of 
the heavy Majorana neutrino density to entropy density ratios Y/v on temperature. 
One can see clearly that for high temperature region z < 0.1, Y/Vi is nearly equal to 
Yjsr 2 - However, as z becomes larger, the asymmetry from the heavier right-handed 
neutrino decays drops quickly, thus, the lighter one gives dominate contribution. 

We plot the evolution of the lepton family asymmetries in Fig. [51 Fig. and 
Fig. El By investigating the family structure of these figures, we observe the follow- 
ing features. 
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(1) The total lepton asymmetry 
Yi, increases with larger z and reach 
its maximum at z ~ 0.6. The 
behavior of Yl in small z region 
can be understood as follows. Since 
lm[(yly v ) 2 ]i2 = -Im[(yjy^) 2 ] 2 i, the 



term ^ 



J 12 



y, 



1 e? 7 £, ? in (|H^5l is pos- 



itive for k = 2 while negative for k = 1. 
In addition, for fl<l, from Eqs. . 
(|3-14|) and (|3-7[) , one can obtain that the 
order of ^ To is the same for k = 1 

i 

and k = 2. However, around z ~ 0.1, 
the deviation from equilibrium distribu- 
tion is larger for the heavier Majorana 
neutrion Nr 2 . Thus, the magnitude of 
the positive one is much larger than the 
negative one. Then, Yl increases in the 
small z region. This feature is in contrast 
to the previous worlPJ where only con- 
tribution from the lighter right-handed 
neutrino are taken into account. As 
z becomes larger, the positive contribu- 



te 




1 10 

z=M/T 



Fig. 5. The lepton asymmetry Yl {Yl) from 
Nr 1 (Nr 2 ) decay. The total lepton asym- 
metry is denoted by Yl. 7l = 0. The 
shaded part shows the bounds on Yl in 
14-6H via experimental measurement^^ 
at 90% C.L. 
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tion to Yl from the heavier Majorana neutrinos decreases and the contribution from 
the lightest heavy Majorana neutrinos decay becomes dominant. We can see that 
the sign of Yl changes at the intermediate regions (1 < z < 10). Then the lepton 
asymmetry from the the heavier Majorana neutrino (Nr 2 ) is compensated by the one 
from the lightest heavy Majorana neutrino (N^). This feature can be seen clearly 
by showing the contribution to lepton asymmetry from N Rl and Nr 2 decays sepa- 
rately, as desplayed in Fig. OD Yl will further decrease and it tends to be a constant 
asymptotically at low temperature. It is important to include the contributions from 
both heavy Majorana neutrinos for the lepton asymmetry and its evolution. 




z=M/T 

Fig. 7. The evolution of the lepton asymmetry and lepton family asymmetries with the CP violating 
phase jl = f • 

(2) In Fig. H3 we have shown the lepton family asymmetries for the two different 
choices of CP violating phase jl- The figure on the left corresponds to 7^ = 
and that on the right corresponds to jl = Despite of similarity in their shapes, 
the dominant contribution to Yl comes from Y^ for jl = 0, whereas, Y T in case of 
"fL = 7r - The results indicate the lepton family asymmetries Y^. are sensitive to the 
CP violating phase jl- To understand about the feature (2), let us focus attention 
again on term in (|3-35j) which is propotental to Im[(y v y v )i2(yu)i\(yu)i2\- Since 
the matrix y\,y u is not related to the matrix Ul, and 7^ just appears in matrix Ul 
(|2-3|) . the 7i-dependent terms in sfj^j only come from the quantity {y v )*i{yy)i2 in 
(f3~T3|t . Using (J23J) and (g3J), we otbain 



{yy)*l{yy)i2 



{{ml 



rrin) sm( 



■ cos 



0r{\(UlU 2 -\(UlU 2 ) 



+m 2 m 3 [cos 2 e R (U L )UUL)i3 - sin 2 R (U L ) a (U L )* a ] } . (4-9) 



Note only the second term of ()4-9|) is relevant to CP phase 72,, and is proportional 
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to {Ul) i2 (U L )is given by 



(U* L )i2(U L ) 13 = i sin 20 lis sin £ ' L12 e^-^ 



(C/l)22(^l)23 



(^) 32 (C/l)33 



1 



sin 20li 3 sin 6» L i2 sin 2 0L23e 



+ cos 6*1,12 cos 6 Liz sin 26*1,23) e i7i ^ 



76 



«7L 



1 



sin 20L13 sin 9 L i2 cos 6> L23 e 



-iS L 



cos 6» L i2 cos 6» L i 3 sin 26*1,23) e 



1<1L 



,;, 



*7l 



(4-10) 



From (|4-9j) and (|4-1()|) that in cases of 7^ = 0, 7r, the 7^-dependent terms of e\ 7^ 
have opposite sign, leading to that is dominant in the total lepton asymmetry in 
case of 7x, = 0, while Y T is dominant in case of 7^ = vr. In Fig. [71 we have shown the 
evolution of the lepton family asymmetries for another choice, 71, = In this case, 
\Y^\ and \Y T \ are much larger than Yl- In contrast to the total lepton asymmetry, 
the lepton family asymmetry, for example, Y^ from Nr 2 decay is dominant and can 
not be upset by the contribution from decay. 



200 
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100 

v so 

I 

■50 
-100 
-150 



-■ *- 



-* ■ 
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y L (degree) 

Fig. 8. The lepton family asymmetries as a function of the CP phase 71,. The data points denoted 
with star are extrapolated. 

The understanding is illustrated in Fig. |Bl We show the lepton asymmetries at 
the low energy and study the dependence of the asymptotic (z 3> 1) values on the 
CP violating phase jl- In contrast to case of Y^ and Y T , Y e and the total lepton 
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asymmetry Yl is insensitive to the value of CP phase jl and is nearly constant. 
In addition, we find from this figure that for the inputs we chose, the size of the 
lepton family asymmetries Y^ and Y T can be of order 10 -8 and cancel each other 
mostly, leaving the total lepton asymmetry of order 10~ 10 which is consistent with 
the experimental observations.^^ 

§5. Conclusions 

In this work we have presented a detailed calculations for primordial lepton 
family asymmetries in seesaw model by taking into account the effect of chemical 
potentials for all SM particles. The minimal CP violating seesaw model is used in 
obtaining numerical results. The main results of our paper are as follows. 

(1) Subject to constraints from neutrino oscillation experiments and observed cos- 
mological baryon to photon ratio, the order of the total lepton asymmetry Yl can 
be different from that for each generations Y^. . There exist scenarios in the minimal 
seesaw model, both K, and Y T are of order 10 -8 , whereas Yl is of order 10 -10 due 
to the cancellation. 

(2) The lepton family asymmetries are sensitive to CP violating phase m the 
matrix Ujj. In contrast to this, the dependence of the total lepton asymmetry on jl 
is quite weak. Note that the total lepton asymmetry is sensitive to the CP violating 
phase 7/? in V/PJ . 

(3) Both of the heavy Majorana neutrinos give important contributions to the pri- 
mordial lepton family asymmetries. The heavier Majorana neutrino decay can dom- 
inate in the lepton family asymmetry. 

The new features we found provide useful informations and may be observable in 
future cosmological experiments. As a future extension of the research, it would be 
interesting to investigate the correlation between CP violation for the lepton family 
asymmetries and the low energy CP violation of neutrino oscillationsP^''^'^'^'^ 
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Appendix A 

Source terms and C nk 



In this Appendix, we show the definitions of the source terms C;. and C nk . 
Writing the source term as CV = CP + CP, CP and C, stand for contributions from 
the decays and the inverse decays of the heavy Majorana neutrinos and two-body 
scatterings, respectively. CjP is defined as, 



cf = j dn 1 dn 2 dn N (27r) 4 5 4 ( Pl + P2 



Pn) 
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E [\M(N Rk (s z )^lr<p+)\ 2 (1 + - /,_) 

k,s z 



+ |A^(iV Rfc ( S ,)^^ )| 2 (1 + /»(1 - /„.) 



|A^(iV Rfc ( S ,)^/+r)| 2 (l + V)(l-/ /+ 



M(N Rk (s z )^V i( f> *) (i + / ^)(l-/ F J 
+ (1 - MO) [/*-/,+ \M(ltr^N Rk (s z ))f 

2 

+ /<AO*/Fi 



-U + f lT \M(l-^N Rk (s z ))\ 2 

-/^/^iMM^iv^^))! 2 ]}, 



(A-l) 



where dilj = ^^Ie- and |A4(a + b + • h + j + • • -)| z stands for the amplitude 

squared for the process a + b + • • — »i + i + • • 

The contributions from two particles scattering processes can be divided into 
four parts, 

(A-2) 



(1)R 



+ d 2)R + cr iL + c; 



(3)fl 



,(4)i? 



The terms c/ x)i? (x = 1, 2, 3, 4) are defined as follows, 



c h )R = J dn 1 dn 2 dn 3 dn^2TT) 4 5 4 ( Pl + P2 -P3- pa) x 

E [MO/^IM^O^0l)| 2 (1 - /,r)(l - ^ 



k.s z -,cx- 



+f N {s z )h l \M{N Rk {s z )bt^lrt a R )\\l - /,_)(1 - / t «) 

-/7v(s,)/ iS |^(iV Rfc ( S ,)^ + ^)| 2 (l " /,+ )(! " /tg) 
-/iv( S2 )/5«|M(^ fc (^)6^^)| 2 (l - /,+ )(l - 

+/iv(^)/rJ^(^(^)^^)| 2 (l - - /s) 



+/Ar(^)/ t «|M(7V Rfc ( s ,)t2^^)|^(l - f Vi ){l - f t% ) 
-fN(s z )f t% \M(N Rk ( Sz )t a R ^t a L )\ 2 (l - / Fi )(l - /ta) 

-/^(^)/^|^(iV Rfc ( S ,)t^I7^)| 2 (l - f Vi )(\ - /ja) 

-4r/^N(0L^ fc (^)^)| 2 (l - - 

-/ir/tsi^tar^^^jfeS)^! - .m^))(i - 

+/z+^|-M(/+^-^(^)^)| 2 (l - fN(s z ))(l - ft,) 

+fi t hi\M{itr R ^N Rk { Sz )b a L )\ 2 {i - f N ( Sz ))(i - f- bt ) 

-fuJr L \M{vr L ^N Rk {s z )T R )\ 2 (i - f N (s z ))(l - ft*) 
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-UJi R \M(vit%-+N Rk (s z )tl)\ 2 {l - f N(Sz) )(l - hi) 
+fvJt 1 \M(v t tl^N Rk (s z )t a R )\\l - f N (s z ))(l ~ ft%) 

+hJt % \M(vr R -^N Rk (s z )r L )\\i - f N (s z ))(i - fr L ) , 

a is color indices of quark. C^ R (x = 2,3) are given as follows, 

c u )R = j dn 1 dn 2 dn 3 dn^2Tr) 4 5\ Pl + P2 - P3 - pa) x 

£ [fN(s z )f lt \M(N Rk (s z )l+ - b a L t a R )\ 2 (l - f t «)(l - f n ) 

k,s z ,a 

+fN{s z )U i \M{N Rk {s z )u i -+ ttt a R )\\l - f- tt ){l - f t% ) 
-fN^frlMiN^l- - btt a R )\\l - f bt ){l - f {% ) 
-fN(s z )UM(N Rk (s z )u t - t a L i a R )\ 2 (l - / t g)(l - fy*) 
-ft % f n \M(btt R -+ N Rk (s z )lt)\ 2 (l - fN(s z ))(l - /,+ ) 
-ftzft % \M(t a R t a L -+ N Rk (s z )i7 i )\ 2 (l - f N (s z ))(l - M) 
+fbtft« R \M(btt R -+ N Rk (s z )l-)\ 2 (l - f N (s z ))(l - /,_) 

+h R ft t \M{tli R - N Rk {s z ) Vi )\ 2 {\ - f N (s z ))(l - fj . 

cf )R = J dn 1 dn 2 dn 3 dn^2Tr) 4 5\ Pl + P2 - P3 - pa) x 

< £ [/^+|-MG+<T-^ + )| ,2 (i + V)(i - %) 

-U + fi 7 \M(lJ<f> + ^lt<f>- )|' 2 (1 + - /,+ ) 

+U-f lt \M(lt<l>- - ^/)| /2 (1 + / o)(l - /„,) 
-fpU^M^ - /+r)l' 2 (l + /*-)(! - / z +) 
+/0O./ Pi |M(^°* - ^/)| /2 (1 + - /„.) 

-fifif^MM -+ ^°*)|' 2 (1 + - /pj 

+f 4 fi-U i \M{p i <t - ^-0 + )|' 2 (l + - %) 

-U + f l7 \M(lJ<l> + - P^°*)|' 2 (l + />,)(1 - /pj 
+U-fi+\M(l^<f>- - C0 + )I' 2 (1 + /*+)(! - /,-) 

-/*+/,- l>W ^ )l' 2 (l + " /,t) 

+U-f lt \M(l}<l>- - ^°)| ,2 (1 + - /„,) 

- /+<nr 2 (i + - /,+) 
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+f^fu J \M(u ] 4P* - ^/)| /2 (1 + j»(l - /„.) 
-Uof Ut \M(^ - ^°*)|' 2 (1 + V)(l - /p,) 

- ^0 + )| /2 (i + U + )(i - f^) 
-f^f^Mil'^ - ^/*)|' 2 (1 + /^)(1 - h 3 ) 

+E - ^-)i ,2 (i + - /,+) 



-U-fMM{i+</>- 

3 J 




k> 


+f<t>-fi+\M(itr 


^^°*)|'2(1 + /^)(1 


-h) 




+ /,_)(! 




+f^fn\M{n<jP* 


^^ *)|' 2 ](i + V)(i 


- fpj) 




^/*)l'2(l+/^)(l 


-M) 


+f^f Vj \M{ujctP -+ l-<f> + )\' 2 (l + U + )(l - 


f l7 ) 




-^°)r 2 (i+/ o)(i- 


fvj) 


+U + fr\M(lJ(f> + 


-/r0+)|' 2 (l + / 0+ )(l- 


-/,-) 


-U+fi-\M(ir4> + 


-/70 + )r 2 (i+/ 0+ )(i- 


-v 


+Uo*M\M(vi<t> * 


— z/0->r 2 (i+/ -)(i- 


-M 


-u-f lt \M(i+r 


^^ o *)i /2 (i+/^)(i 


-M) 



+f 4>0 f Uj \M(v J 4>° - ^°)r(i + /»(i - /,j 

SMMiyrf -+ ^°)| ,2 (1 + /»(1 - /„,) 
+U + f l7 \M(lJ<f> + - ^°)| ,2 (1 + j»(l - /„J 

- /70 + )| /2 (l + V)(l - /,-)] } , (A-5) 

where, 

\M(a + b^c + d)f = \M(a + b ^c + d)\ 2 - \M(a + b^N Rk ^c + d)\ 2 . (A-6) 
Finally, C^ R is given as, 

)fl = y dn 1 dn 2 dn 3 dn 4 (27r) 4 5 A ( Pl + P2 -p 3 - pa) x 

< ^ E [ - ^)i 2 (i - /,j(i - /„.) 

/,../,.>(•>.'', - /*0°*)| 2 (l + />,)(l + V) 

+u-u-\m(<i>-<i>- - i;ij)\ 2 (i - /, r )(i - %) 
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-fir fr \M(l~l- - <T<m 2 (l + /*-)(! + V) 
+/pJp j |M(^ - W)| 2 (l + /^o)(l + V) 

- - ^)| 2 (i - - / p .) 

+/*+/ ;+ IM^ + - <^ + )| 2 (l + /*+)(! + V) 

-u+faiM^^ - /+/;)i 2 (i - /,+)(i - /,+/ 
+E [/*-/«h - ^-)i 2 (i - //-)(! - Ui) 

i 

I, /.-. Mil, v } -+ <T/*)| 2 (1 + U-)(l + fa) 
+fyfv j \M{ltv j - ^V)| 2 (l + + /*°) 

-fafaM{<t> + <t>° - /+^)| 2 (1 - - f Pj ) 

+U-fa\M(<f>-<i> * - ^)| 2 (i - /, 7 )(i - /„<) 

-hjUMilJvi - 0-0°*)| 2 (l + /«-)(! + V) 
+f lt M\M(lp t -+ 0+^°)| 2 (l + / 0+ )(l + j» 

- fafaM(<i> + <f> - Z+^)| 2 (l - - 

+ E [WIM0 + <r - C^)l 2 (! - /r)(! " /it) 

/, // >';/, // - </> + <ni 2 (i + V)(i + fa) 

1 3 

+fitfi 7 \M(lp- - <^ + 0-))| 2 (l + fa)(l + fa) 

-fafa\M(^ + r - ip;)\ 2 (i - fa(i - fa 

+U°fa\M(4>°<t> * - ^)| 2 (i - - /p,) 

-UMMiyiVj - <AV°*)| 2 (1 + / o)(l + fa) 
+U j faM{v j n - 0>°*)| 2 (1 + / o)(l + / o,) 
-fyfalMtt ^ - ^)| 2 (1 - - A) 
+fafa\M((f> + <l> * - ^ + )| 2 (1 - fa(l - fa 

-fa\+\Minl+ - </>V*)| 2 U + fa){l + fa) 
+fa\- \Mivil- -+ <T/)| 2 (l + fail + fa) 
-Uo fa\Mi4>°4>- -^7)l 2 (i - - %) 

+/«<>/,- - ^Dl 2 (i - /, r )(i - 4) 
-/, r / Pj .|M(^ir ^ 0V")l 2 (i + /*-)(l + V) 
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-/^/■H-^V* - )l 2 (i - - /,+)] } • 



(A-7) 



In deriving (|A-7|) . the statistical factors are counted as follows. For the first 
eight processes with i = j, we need to multiply a factor 2 compared with those i ^ j 
since i = j cases corresponds to |ALj| = 2 processes. However, there is a symmetric 
factor | in the case of i = j, whereas \ f° r other cases. Therefore we have a common 
statistical factor for i = j and i ^ j. 

Next we calculate C Uk in which can also divided into decay part C® k and 
2^2 reaction part C^ k . Similar to the calculations of the C\ t , 



eg = J dn 1 dn 2 dn N (27T) 4 5\p 1 +p 2 -PN) x 

£{(!-/*(**)) [\M(if<p- ^N Rk { Sz ))\ 2 U-f lt 



+\M(Pi<t>°* - N Rk ( Sz ))\ 2 f^f Pi + -> AT flfc ( s ,))| 2 /^ + / ir 

+ \M(u i( t>°^N Rk (s z ))\ 2 f^f Ut ] 
-f N (s z ) [\M(N Rk (s z ) -> ^°)| 2 (1 + - /„,) 

+|M(i\fe fc ( a ,) - Zr0+)|2 ( i + ^ +) (i _ /{ _) 

+ \M(N Rk (s z ) - Z+<T )| 2 (1 + - /,+ ) 

+ 1^^^) - z^°*)| 2 (l + - /pj] } • (A-8) 

Since the smallness of Yukawa couplings of lepton and light quark, only contributions 
from the processes involving top (anti-top) quark to the reaction part are important, 
i.e., 



C, 



R 



>i k 



(A-9) 



with, 



c {i)R 



J dn 1 dn 2 dn 3 dn 4 (2ir) A 5 A (p 1 +p 2 -p 3 - pa) x 

£ [firh t \M{lfbt -» A fifc (^)tl)| 2 (l - - /%) 



+/ir/tsl-W(T*S 

+f lt fb* L \M(lfbt 

+f lt ft % \M(l+i a R 

+fuJt R \M(t a R ^ - 

+h i h< k \M{t a B p i - 

+fuJtt\M(t a L iy t - 

+fuJt 1 \M(t a L v l - 



,N Rk (s z )bt)\\l-f N (s z ))(l-f bl ) 
,N Rk (s z )t a R )\ 2 (l-f N (s z ))(l-f tR ) 
, N Rk {s z )bt)\ 2 (l- f N {s z )){\- f n ) 
-N Rk {s z )tt)\ 2 {l-f N (s z )){l-f tl ) 
N Rk (s z )it)\ 2 (l-f N (s z ))(l-f- tl ) 
N Rk {s z )t a R )\ 2 {l-f N {s z )){l-f- taR ) 
N Rk (s z )t a R )\ 2 {\- f N (s z ))(l- f t% ) 
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and, 



-fN(s z )ft t \M(N Rk (s z )i a L -+ l-b a R )\\l - - f n ) 

-f N {s z )h t \M{N Rk (s z )bt -+ ltm\ 2 (l ~ /,+ )(! " ft" R ) 
-fN{s z )fti\M{N Rk {s z )tt -+ t a R n)\\l - ~ ft % ) 

-fN(s z )ft t \M(N Rk (s z )it -+ ^)| 2 (1 " /pJ(1 " fi % ) 
-fN{s z )f- tl \M{N Rk (s z )i a R -+ itn)\\l - f Vi ){i - hi) 
-fN(s z )ft«\M(N Rk (s z )t a R -+ ^)| 2 (1 " /pJ(1 " ftt) 

-fN{s z )h t \M{N Rk ( Sz )bi - i m -*s)i 2 (i - Ar)(i - /&) 

— fN(s z ) ft<^\M.(N Rk (s z )t R - Z+62)| 2 (l - /,+ )(! - 



= y dn 1 dn 2 dn 3 dn 4 (2ir) 4 5\p 1 +p 2 -p 3 - P4 ) x 
E ^ iv Rfc (^)/+)| 2 (i - - /,+) 

\M{bffi - iV Rfc ( Sz )/-)|2 ( i _ /7V ( Sz ))(i - /,_) 
+/f«/tgl^(^ - iV Rfc ( Sz )^)| 2 (l - f N (s z ))(l - f Vi ) 
+ftift % \M(t a R it -+ N Rk {s z )n)\ 2 {l - f N (s z ))(l - f Pt ) 
-f N (s z )f lt \M(N Rk (s z )l+ -+ b a L t a R )\ 2 (l - / t « )(1 - f n ) 
-fN(s z )f lT \M{N Rh (s z )lT - b a L i%)\ 2 (l - f bl )(l - f%) 
-f N (s z )UM(N Rk ( Sz )u t - t a L i%)\ 2 (l - f lt ){l - f t% ) 

-fN{8 s )fv i \M{N Rk {8 z )v i - ^)| 2 (1 - /f«)(l - /*a) 

Appendix B 

Reduced cross sections a(s) 



'•Kl 



a 



N,m 



1 

2vr 



(A-10) 



(A-ll) 



The reduced cross sections with N Rk exchange can be expressed as, 



(B-l) 

where m denote the types of the contribution and take m = 1, 2, 3, s, il, t2 respec- 
tively, ij denote the types of lepton family. After straightforward calculations, we 
obtain, 



sak 



2(s - a k ) 2 + c k a k 



2ak_ 



Ok 

s 



1 + 2 



s + a k 



In 1 + 



(B-2) 
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+ 



(L + l) 



s 1 1 

2DiD k + D~i + D k ~ 
1 1 



+ 1 + 



Ok 



&k — 0,1 D k 
1 



In 1 + 



ai -a k DiJ 



- -M h | I + 



a; 
s 

o.k 



fkk 

In, 2 



flk 
JN,2 



fkk 
JN,3 



flk 
JN,3 



fkk 
JN,s 



fkk 

JN,n 



fkk 
JN,t2 



flk 
JN,t2 



2s 



+ 



4a k 



In 1 + 



s + a k I + 2a k 
(a; - a fc )(s + a; + a fe ) 



OA: 



(s + 2a/)ln( 1 + — 



-(s + 2a fc )ln 1 + — 
ai 



1 



sa fc 



2 (5 - a k ) 2 + c k a k 



\/aia k 



1 



s 2 



2 (s - a fc ) 2 + c k a k 
s 



Ak 
JN,i 



D,D 



Ak _ 2 ^QiQk 
s 

O'k 



s + a k 

- 2 +( 1+2 !W 1 + £)' 



i^k 



In 1 + 



a. 



-In 1 + 



-2 + 



az - afc 



«,(l. + ^)ln( l+- 



(l + ^)ln(l + A 



The corresponding reduced cross sections are denned by, 



a 



N,3 



(7 



(B-3) 
(B-4) 

(B-5) 
(B-6) 
(B-7) 
(B-8) 
(B-9) 
(B-10) 

(B-ll) 

(B-12) 
(B-13) 
(B-14) 



For the s channel Higgs exchange process — >t b a , the reduced cross section 

reads, 



*fc ^^Hs)(N R J- - = ^ (y\y t ) \(y u ) lk f (l - , 



(B-15) 



where iV c is the color number of quark , y t is top quark Yukawa coupling , and for 
the t channel Higgs exchange process Nn k t a ^l^b a , 

aft = 4j2Hm Rk t a - a ) = § (ylvt) \(vM 2 (i - f ) 
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s — 2d k + 2clh (ik — 2a# an 

in ■ 



S - dk + OH S - dk S -dk+CLH 



(B-16) 



where a&, s are denned in Q3-4j) and (|3-5|) . respectively, a// = Mfj/Mf with m# 
being the Higgs mass, c k = (T£/Mi) 2 , = 4^r fei , 

i 

1 S — dh . 

— = jz ^r 1 , (B-17 

t>k {S-d k ) Z +d k C k 

is the off-shell propagator of Nji k . 
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